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We analyze the neutrino induced charged current coherent pion production at the energies of 
interest for recent experiments like K2K and MiniBooNE. Medium effects in the production mecha- 
nism and the distortion of the pion wave function, obtained solving the Klein Gordon equation with 
a microscopic optical potential, are included in the calculation. We find a strong reduction of the 
cross section due to these effects and also substantial modifications of the energy distributions of 
the final lepton and pion. 
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Recent experimental results from the K2K collaboration show a significant deficit of muons in the forward scattering 
events with respect to the simulations, which could limit the accuracy of the predicted neutrino energy spectrum at 
the far detector This deficit could be due, among other possibilities, to the overestimation of charged current 
(CC) coherent pion production. A later search Q found no evidence of CC coherent pion production in 12 C and 
\q • obtained an upper limit for the fraction of this process to the total CC interaction well below some estimations based 
£NJ ] on the Rein and Sehgal model of Ref. [3|. Preliminary MiniBooNE results for neutral current 7r° production show 
a similar deficit [H in the forward direction when compared with different MonteCarlo (MC) models: NUANCE [f|, 
NEUGEN @ and NEUT [3] . The discrepancies between these MC simulations indicate the considerable uncertainties 
in the theoretical description of coherent pion production. In addition, MiniBooNE has collected a large set of data 
for 7r + production in 12 C induced by muon neutrinos of energies around 0.7 GeV p|. A fraction of these pions is 
created coherently, so that a realistic description of the coherent process is required to analyze and understand these 
_, ■ data. 

CC coherent pion production has been observed experimentally at higher energies and for several nuclei [t| [l(J HH, 
Eli EE G3- These experiments were studied theoretically using models based on PCAC [1, 0j| which described well 
the results. Work along the same lines has been carried out at low energies in Ref. [l^]. The only nuclear medium 
effect considered in these calculations is the distortion of the final pion. 

Other approaches have tried to incorporate some additional nuclear medium effects that modify the weak pion 
production. In Ref. [l7T |. the authors used the impulse approximation with undistorted pion waves, but already 
modified the A resonance properties, and therefore the production mechanisms, using an effective A mass. The 
importance of these nuclear effects was demonstrated there by comparing the results with those obtained using a free 
A. Kelkar et al. [HI developed a more sophisticated treatment of the A in the nuclear medium, and also included the 
, ' final pion distortion by solving the Klein Gordon (KG) equation with a pion nucleus optical potential. This model 
predicted a very low cross section, compatible with the recent results of Ref. 0. Nonetheless, there were several 
approximations in this work. On the one side a non relativistic reduction of the hadronic current was done, on the 
other side all the transverse parts of the amplitude were neglected. Whereas these approximations are quite reasonable 
to get an estimate of the cross section for this process, as discussed in Ref. 1181, a more complete calculation is required 
now that new data are becoming available. More recently, Singh et al. \19L l20( used similar medium effects on the 
production mechanisms and improved on the description of the elementary v + N — * N + fjT + 7r + process by using 
a fully relativistic calculation of this process and including all pieces of the amplitude. However, the pion distortion 
was implemented in the eikonal approximation, which is not very reliable for the low energy pions that apparently 
dominate this reaction. Our aim in this paper is to improve the model of Kelkar et al. [l8| in a similar manner, 
using a more complete and relativistic elementary amplitude but still keeping a more realistic treatment of the pion 
distortion which is calculated solving the KG equation. 

In the following, we describe the formalism, including A production and decay, medium effects on the production 
mechanism and the optical potential responsible for the distortion of the pion. In Section [ill! we present our results, 
compare with the available experimental data and make predictions for other nuclei and observables. Finally, our 
summary and conclusions appear in Section [TV] 
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II. THEORETICAL MODEL 



In the CC coherent pion production induced by neutrinos (y + A — ► A + A* - + tt + ), the nucleus remains in its 
ground state. The process consists of a weak pion production followed by the strong distortion of the pion in its 
way out of the nucleus. First, we must consider the elementary process of 7r + production (y + N — > pT + N + n + ). 
There are large experimental uncertainties in this cross section at the low energies relevant to this work, which can 
be clearly appreciated by comparing, for instance, Refs. [2l[ and (22J. These discrepancies are ultimately translated 
into different values for the A resonance axial form factors that appear in the theoretical models, and are a source of 
uncertainty for the calculation of the pion production cross sections. In any case, it seems to be clear that for neutrino 
energies below 2 GeV. this process is dominated by a A(1232) excitation, v + p — ► pT + A ++ or v + n — > p~~ + A + 
followed by its decay [23|, [24J, [25|, ]2Q, 127], l28|, [29| , and therefore it is larger for protons than for neutrons as discussed 
below. The situation is different for ir° production or at higher energies. In that case, the A mechanism is not enough 
to provide a good description of data [291 ] . After the introduction of the elementary model for the reaction we discuss 
how it is modified in the nuclear medium, due to the density dependent changes of the A resonance peak's position 
and width. Finally, we study the distortion of the pion wave function using a pion nucleus optical potential. 

A. A production and decay 

The matrix element for the elementary process v^{k) + n(p) — > A + (p') + p~(k') is written as 



with the leptonic current 



and the hadronic current 



G 

M n . A + = cosd c l a J£ , (1) 



l a =«<(fc')7a(l -7s)u,, t (k), (2) 



J&=W)A" a u(p) ) (3) 



H^rig^h - 9 M 7 Q ) + %(g» a q ■ p' - «V°) + c 5 V Q + , (4) 

where M is the nucleon mass, ip^p') and u(p) are the Rarita Schwinger and Dirac spinors for the A and the nucleon 
of momentum p' and p, q — p' — p — k — k' is the momentum transfer, Cj and Cf (i = 3, 4, 5, 6) are the vector and 
axial vector transition form factors. The amplitude for the process on protons, ty^fc) + p(p) — > A ++ (p') + /i~(fc'), is 
related to the previous one by an isospin factor 

M pA+ + = ViM nA + . (5) 

The vector form factors can be related to the electromagnetic ones. The conservation of the vector part of the 
current implies that C§ = 0. The assumption of M\ + dominance for the A electroproduction amplitude gives [23[ 



and for we take [3C 



c v = 2 - 05 (7) 

(l-g 2 /0.54GcV 2 ) 2 ' 



Except for C§ that can be related to C§ using PCAC, there are no other constraints for the axial form factors. 
We use the following ones, fitted to neutrino scattering data [HI [H, lEl S3, [H, H [H H E3| 



Cl 3A Jq 2 )=Cf(0) 



h - q 2 



(8) 
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and 

<W) = <#-^, (9) 

with C^(0) = 0, C£(0) = -0.3, C^(0) = 1.2, a 4 = a 5 = -1.21, 6 4 = &5 = 2 GeV 2 and for the axial mass we take 
A/a = 1.28 GeV [32]. This set of form factors produces a good agreement with weak A production data induced by 
neutrinos on nucleons [26j . 

A new analysis of world electron scattering data [38| allows to go beyond the M\ + approximation and update the 
vector form factors [39(. This information calls for a new extraction of the axial form factors. Some steps in this 
direction have been recently taken in Refs. [27],[4(|, using models for the elementary reaction which include, apart 
from A excitation, some background terms. For the sake of consistency with our description of pion production on 
the nucleon explained previously, based on the A dominance, we stick to the set of form factors given above. 

We still need to take into account the A decay into a pion and a nucleon. We use the following Lagrangian to 
describe the ANn transition [4l| 

Ca*n = — *pTt • * + /i.e. , (10) 

where is the isospin 1/2 to 3/2 transition operator and 4> and \& are the A, pion and nucleon fields respectively. 
The isospin operator produces a factor 1 for the decay A ++ — * 7r + + p and a factor for the decay A + — > tt + + n. 

This, together with Eq. (5), implies that the amplitude for CC ir + production on the proton is three times larger than 
on the neutron. The coupling constant /* = 2.13 is such that the experimental A — > ttN width is reproduced. At the 
high A invariant masses reached with the K2K or MiniBooNE neutrino energies, the finite size of the hadrons becomes 
relevant. This can be empirically taken into account with a form factor F(p') that modifies the ANn coupling. Here, 
we adopt the following ansatz 

p w = A 4+( /-Mir 

with A = 1 GeV, used in coupled channel studies of pion and photoproduction of baryonic resonances The new 
hadronic current which should replace J% in Eq. (JTJ), already incorporating the decay into a pion and a nucleon, is 
given (for the case of A + ) by 

JL = -^=£- P :F(p')u(p f )D(p')A a0 A^u(p) , (12) 
v 3 Tn n 

where p^ and pf are the pion and final nucleon momenta so that p' = p v + pf. The A propagator is given by 

D{pl) = (W + M A )(W -M A +iT A /2) ' (13) 
where W = \fp' 2 - The energy dependent A width is 

1 ^/'V^^AJ 



r A = 7T — Fip'f^plcm, (14) 

with Pn^cm the pion momentum in the A rest frame. Finally, the spin 3/2 projection operator is given by 



= - (V + Ma) [ Bafl --— + - 7q7/3 ) . (15) 



B. A in the nuclear medium 



The A properties are strongly modified inside the nuclear medium and have been the subject of intensive study, 
both experimental and theoretical, for many years, see i.e. (43l l45l l4rl WA l48l. l49l. [50j and references therein. In this 
work we use the results from Refs. [50l[5l| where the A selfenergy is calculated in a many body approach as a function 
of the local baryon density p(r). This model has been extensively tested in pion induced inclusive processes [52l [53|. 
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photonuclear reactions [Hil . l55l l , e lectron scattering [56| and even for coherent pion production induced by photons, 
electrons or nuclei [57l. l58l |5& l60l]. In the nuclear medium the A resonance acquires a selfenergy because of several 
effects such as Pauli blocking of the final nucleon and absorption processes: AiV — > NN, AN — > NNtt or ANN — > 
NNN. The real part can be parametrized as 



Re£ A (p)=Re£i(p) + ^ (£- 
9 V m.,r 



9 P ' 



40 MeV — , 
Po 



(16) 



where po — 0.17fm~ 3 is the normal nuclear density. Here, in addition to the attractive proper selfenergy ReS^(p), 
the effective repulsive contribution, that comes from the iterated A-hole excitation driven by the Landau Migdal 
interaction with g' = 0.63, has been added [H3 |. The imaginary part is parametrized by the expression 



ImE A (p) = C Q ^. 

Po 



CA2 ( — 

Po 



CA3 ( — 

PO 



(17) 



where the terms with the coefficients Cq, Cai and Caz correspond to the processes AN — ► NNtt, AN — ► NN and 
ANN — > NNN respectively. The values of Cq, Cai, Ca3, a, (5 and 7 can be found in Eq. (4.5) and Table 2 of 
Ref. (52| . The parameterizations are given as a function of the kinetic energy in the laboratory system of a pion that 
would excite a A with the corresponding invariant mass, and are valid in the ran ge 8 5 MeV < T v < 315 MeV. Below 
85 MeV the contributions from Cq and Ca3 are rather small and are taken from [5l|, where the model was extended 
to low energies. The term with Cai shows a very mild energy dependence and we still use the parameterization from 
Ref. [52| even at low energies. For T n above 315 MeV we have kept these selfenergy terms constant and equal to 
their values at the bound. The uncertainties in these pieces are not very relevant there because the A — > Ntt decay 
becomes very large and absolutely dominant. Finally, the Pauli blocking of the irN decay reduces the r A free width 
which now reads as 



r5 



(18) 



The angular integrals I\ and I2 can be found in Appendix B of Ref. [5l| . This selfenergy is taken into account by 
making the substitutions M A — > M A + ReS A and T A /2 — > r^ auli /2 - ImS A in the propagator of Eq. JTc 



C. Cross section 



The hadronic current is further modified in the nucleus, where the nucleons are bound and thus have a momentum 
distribution. In the impulse approximation, and after summing over all nucleons, we can write it as 



2 



1 / ifireKG-P*)-? 



, s . Pn(r) 
Pp( r ) + ~^r~ 



V3^F(p')D(p')p«Tr{u(0)A af3 A^u(0)} , (19) 



where D(p') is the in-medium A propagator and the trace corresponds to the sum over the nucleons' spin. In the 
evaluation of the trace we have taken an average nucleon momentum for the spinors neglecting corrections of order 
(p/M) 2 . Given the large mass of the nucleus, the pion energy is taken to be equal to q° = E v — and this defines 
the modulus of the asymptotic pion momentum. On the other hand, a momentum of q — is transferred to the 
nucleus. In the evaluation of the amplitude we assume that this momentum is equally shared by the initial and the 
final nucleon, so that they have pi = (pV — q)/2 and pf = (q — p^)/2 respectively. This prescription has also been 
used in Refs. [13, HH for coherent 7r° photo- and electroproduction. The approximation is based on the fact that, for 
Gaussian nuclear wave functions, it leads to an exact treatment of the terms linear in momentum of the elementary 
amplitude and allows for a consistent description of the pion- nucleon and pion- nucleus kinematics [61( . Pion distortion 
is taken into account by the replacements 

e-*** - (20) 

and 

f^e-^-r iV^ ut (p;,r), (21) 

where 4>* ut (p^ , r) is an outgoing solution of the KG equation for a pion of asymptotic momentum p^ calculated with 
the optical potential described below. To test the validity of the eikonal approximation we have also calculated the 
cross section with the pion wave function 



££ fc (&,0 = e- ,p -" r e J ' ~^^ dz , (22) 
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where r = (b, z) and II is the pion selfenergy described in the next section. 

The cross section for the coherent process v + A — > A + pT + ir + is then given by 



da 1 |fc'|K| 



dfledEedQ^ 8 |fc| (2tt) e 

with 



\M\\ (23) 



M = cos c l a J^. (24) 

D. Pion optical potential and distorted wave function 

The pion wave function is the solution of the KG equation with an optical potential. Since most of the produced 
pions lie in the energy region around the A excitation the A-hole model can be used. In this model, which we take 
from Ref. [50(, the optical potential is given by 

V opt = £ , (25) 

where w is the pion energy in the laboratory system and the selfenergy II reads as 

M 2 „ V 

II = -47r q 2 - — . (26) 

Here, s is the Mandelstam variable of the irN system and q is the pion laboratory momentum. Finally, V is given by 

_J_/ j T\ 2 f Pp + p n /3 Pn + Pp/3 1 

6tt \m„J \ y/s - M A - Re£ A + ir A auli /2 - i ImS A -^i - M A + 2M - ReS A / ' 1 ' 

where the two terms correspond to direct and crossed A-hole excitations. The neutron densities p n are taken from 
Ref. [62j and the proton densities p p are obtained from the parameterizations compiled in Ref. 63] , both deconvoluted 
to take into account the finite size of the nucleons as done in [62} . The real and imaginary part of the A selfenergy 
Sa and the Pauli corrected width, r A , have been described in section Hi Bl 
In coordinate space and for finite nuclei, this p-wave potential can be cast as 



f2 



M 

2u;V opt (r)=4TT 



7>(r) ^ 1 lj A V(r) 



V - -—A: 



1 + 4Trg'V(r) 2M l + Airg'Vir) 



(28) 



where the first term has the standard Kisslingcr form and the second one accounts for the angular transformation 
from center-of-mass to laboratory variables. The r dependence in V appears via the local density approximation 
p — > p(r). With this potential we solve the KG equation and obtain the pion scattering wave function. The procedure 
is described in detail in Refs. [sH • To asses the quality of the potential one can compare its results for differential 
cross sections with pion nucleus elastic scattering data. We obtain an overall good agreement from light to h eavy 
nuclei at the energies relevant for this work. We can also refer the reader to the article by Garcia Recio et al. [65| . 
where a very similar potential was considered, using as here the local density approximation and the same values for 
the A selfenergy. The main differences with that work are their inclusion of a phenomenological s-wave selfenergy, and 
our inclusion of the A crossed term. We have checked that both pieces produce only minor effects in the A resonance 
region. 



III. RESULTS 



We show in Fig. [T] the differential cross section da/dp^ for CC coherent pion production in 12 C at a neutrino 
energy E v = 1 GeV. This plot shows the effect of the different ingredients of the calculation. The modification of 
the elementary reaction mechanism, through the inclusion of the A selfenergy in the propagator, already produces a 
strong reduction of the cross section. This reduction of around 35% agrees with the results of Ref. [13]. The pion 
distortion further decreases the cross section and moves the peak to lower energies. This reflects the presence of a 
strongly absorptive part in the optical potential around the A resonance peak. The final result shows that the pion 
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FIG. 1: (Color online) Momentum distribution of the coherent pions. 



spectrum is peaked at much lower energies or, equivalently, that the muon energy distribution is peaked at higher 
energies than for the impulse approximation. As expected, the eikonal approximation fails for low and intermediate 
energies, where a better treatment of the pion wave function is clearly required. Our eikonal result differs from the 
one obtained in Ref. [l9| because of their use of the asymptotic momentum in the amplitude, whereas we take the 
gradient of the distorted pion wave function, and our more complete treatment of the optical potential (see Eq. ([26]) 
vs. Eq. (6) of Ref. [H). 

In Fig. [2J we present the same observable averaged over the K2K [6(| and the MiniBooNE (6?| spectra compared 
with the results for a fixed neutrino energy. The consideration of high energy neutrinos widens the pion momentum 
distribution, due to phase space. Even when both spectra are quite different and K2K has a larger average neutrino 
energy (1.3 GeV) than MiniBooNE (0.75 GeV), the peak position stays at the same pion momentum, below the A 
resonance. Furthermore, most of the pions have relatively low energies such that the use of the A-hole model is 
appropriate. 

Although the total cross section and the energy distribution are strongly modified by the nuclear medium, the 
angular distribution of the muons remains relatively unaffected, as can be seen on the left panel of Fig. [3] where 
we compare the result of our model with the impulse approximation, rescaled to match the full calculation at zero 
degrees. On the right panel, the muon angular distributions averaged over the K2K and the MiniBooNE neutrino 
spectra are shown together with the one obtained for 1 GeV neutrinos. In both cases, the consideration of higher 
energy neutrinos leads to a narrower angular distribution. As a consequence, the MiniBooNE angular distribution is 
appreciably more forward peaked than the one at 1 GeV even when its average energy is lower. We have also checked 
that the lower bound on the muon energy used by K2K in Ref. [2] (p M > 450 MeV/c) does not modify appreciably 
this observable. 

We show the total cross section as a function of the neutrino energy and for several nuclei in Fig. 2J with the 
caveat that the pion production model is less satisfactory at high energies. One reason is that mechanisms, other 
than the excitation of the A resonance, become relevant. Also the pion distortion, based on the A-hole model, is not 
appropriate for the high energy pions that can be produced by neutrinos with E v > 2 GeV. Nonetheless, the irN 
interaction is much weaker at high energies than at the A peak and thus distortion effects should be smaller there. 

In Fig. we present the dependence of the total cross section on the atomic number. In this process, the amplitude 
is the coherent sum of the contributions of all participant nucleons. Taking into account the isospin factors, this implies 
that the amplitude is proportional to an effective number of participants defined as P=Z+N/3; here Z and N are the 
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FIG. 2: (Color online) Momentum distribution of the coherent pions averaged over the K2K and MiniBooNE fluxes; the same 
distribution at E v = 1 GeV is also plotted for reference. 

number of protons and neutrons respectively. This could suggest a quadratic dependence of the cross section on P, 
so that for heavier nuclei, the process could be comparatively larger with respect to incoherent tt production or other 
processes. However, there are several reasons that quench the P dependence. First, pion absorption is quite strong 
and forces the reaction to be peripheral. Notice that the cross sections for the inclusive processes are not affected 
by this. Second, the nuclear form factor is narrower for heavy nuclei, and reduces more the contribution from high 
momentum transfers. 

Our result for the integrated cross section for coherent pion production on 12 C averaged over the K2K flux is 
a< coi = x 10~ 40 cm 2 . This value is above the upper limit of 7.7x 10 -40 cm 2 obtained using the ratio between coherent 
and a cc , the total CC cross section, from the K2K collaboration 0] and the value for a cc of their MC calculation. 
Without the experimental threshold for the muon momentum, p M > 450 MeV/c, we obtain — 12. x 10~ 40 cm 2 . 

There are several factors that might conspire to produce the disagreement between our calculation and the exper- 
imental upper bound. First of all, due to the considerable uncertainty in the experimental data for pion production 
cross sections on the nucleon, the axial N-A form factors are not sufficiently constrained. A more complete theoretical 
description of the elementary amplitude with the inclusion of background terms [HI, H(| and heavier resonances [23|, [39| 
could help, but in order to put such a model on a firm ground, more precise data are required. The optical potential 
employed in our calculation is realistic around the A peak (i.e. for pions with 150-450 MeV/c momenta), where most 
of the strength of the reaction actually concentrates, but has room for improvement both at lower and higher energies. 

It is also important to recall that at forward angles, where the coherent process is sizable, nuclear effects play 
an important role and this may affect the experimental separation of the coherent events from the incoherent ones 
which, to a large extent, relies on the theoretical models built in the MC simulations. The situation is illustrated 
in Fig. [5] where we plot the muon angular distributions averaged over the K2K flux (and with > 450 MeV/c) for 
coherent tt + production, together with the main contributions to the total inclusive CC cross section: quasielastic 
scattering (QE) and incoherent A excitation. The calculations of the A part is performed with the same elementary 
amplitudes and in-medium effects as given above for the coherent reaction. For the quasielastic process, we have 
adopted the model of Ref. (68| but updating the nucleon form factors according to [6!| . Nuclear effects include Fermi 
motion, Pauli blocking with a local Fermi gas and the renormalization of the weak transition, which is treated as an 
RPA resummation of particle-hole and A-hole states. These nuclear correlations cause a considerable reduction of 
strength at low q 2 (forward angles), as can be seen in Fig. [51 while they are negligible for cos# M < 0.8. Therefore, if a 
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FIG. 3: (Color online) Muon angular distribution in CC coherent pion production. Left panel: full model vs. impulse 
approximation. The latter has been rescaled to match the full calculation at zero degrees. Right panel: full model at E v — 1 
GeV and averages over K2K and MiniBooNE fluxes. 

model that lacks these correlations is used to extrapolate the data from the region of cos# M < 0.8 to forward angles, 
one might overestimate the QE part, causing an underestimation of the contribution of other mechanisms, like the 
coherent pion production, to the cross section. 

IV. SUMMARY AND CONCLUSIONS 

We have studied CC coherent pion production induced by muon neutrinos + A — > \x~ + 7r + + A. Our model 
takes into account the modification of the production mechanism due to the renormalization of the A properties in 
the nuclear medium and the distortion of the final pion. The distorted pion wave function is obtained by solving the 
Klein Gordon equation with an optical potential based on the A-hole model. Both effects produce a large reduction 
of the cross section with respect to the impulse approximation. The distortion of the pion shifts the peak of the pion 
energy distributions towards low energies. The angular distributions are slightly widened by the nuclear medium 
effects. 

The A dependence of the coherent process has been investigated. We have found that the integrated cross section 
grows more slowly than one would naively expect due to the strong pion absorption and the effect of nuclear form 
factors. Therefore, we do not expect that coherent pion production becomes more relevant with respect to the 
incoherent processes for heavier nuclei. 

We have also studied the cross sections averaged over the K2K and MiniBooNE spectra. In the case of K2K, we 
find a cross section 30% larger than the upper limit estimated in the experiment. One should however remember that 
such upper limit is based on values for the CC cross sections on nuclei (quasielastic, pion production, etc.) which are 
not well known at the rather low energies discussed here. Further improvements in the theoretical description of this 
reaction will require more precise data for the neutrino nucleon cross sections at low energies that could constrain the 
axial N-A form factors. As for MiniBooNE, we give predictions for energy and angular distributions which can be 
useful to compare with and analyze their data. 
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FIG. 4: (Color online) Total cross section for CC coherent pion production as a function of the neutrino energy for several 
nuclei. 




FIG. 5: (Color online) Total cross section for CC coherent pion production at E v = 1 GeV as a function of the effective number 
of participants P=Z+N/3, as explained in the text. 



10 




FIG. 6: (Color online) Muon angular distribution averaged over the K2K flux for incoherent CC quasielastic (QE) scattering, 
A excitation and coherent tt + production. 
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